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properties,  and  that  no  simpler  differential  law  can  describe 
the  observed  features.  In  the  latter  class,  which  defines  a 
viscoelastic  solid  model,  and  possible  initial  or  induced 
anisotropy,  it  is  again  shown  that  no  simpler  differential  law 
can  describe  the  observed  features,  but  that  the  two  types  of 
response  are  not  sufficient  to  determine  the  model.  A  small 
strain  approximation  is  deduced. 
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SUMMARY 


A  frame  indifferent  tensor  relation  between  stress, 
stress-rate,  strain-rate,  and  strain-acceleration  is  constructed 
to  describe  the  response  of  ice  to  maintained  constant  stress 
and  to  maintained  constant  strain-rate  in  uni-axial  stress. 

It  is  shown  that  the  two  types  of  response  reflect  some  common 
properties  of  this  viscoelastic  fluid  law,  but  that  each  type 
contains  further  independent  features,  so  that  both  are 
necessary  to  determine  the  required  response  functions.  The 
shape  of  the  tensor  relation  cannot  be  determined  by  uni-axial 
response.  A  reduced  model  with  response  functions  of  single 
argument  is  constructed,  and  correlated  with  an  idealised 
family  of  constant  stress  responses.  The  extent  to  which  the 
reduced  model  can  reproduce  the  complete  response  is  illustrated 
for  two  different  families.  Next  a  viscoelastic  solid  law 
relating  stress,  stress-rate,  strain  and  strain-rate  is  shown 
to  be  compatible  with  the  same  uni-axial  responses,  but 
dependence  on  the  reference  configuration  through  strain  allows 
a  description  of  initial  or  induced  anisotropy.  A  small  strain 
approximation  is  constructed.  However,  this  minimial  form 
involves  three  independent  response  functions  for  an  initially 
isotropic  solid,  and  constant  load  and  constant  displacement 
rate  responses  provide  only  two  relations.  A  further  restriction 
must  be  imposed  to  complete  the  description^ x^it  the  simplification 
to  constant  modulus  at  constant  strain-rate  is  shown  to  be 
incompatible  with  expected  response.  We  are  now  investigating 
alternative  simplifications.  Note  the  contrast  with  the  minimal 
fluid  model  which  is  overdetermined  by  independent  constant  stress 
and  constant  strain-rate  responses. 
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Viscoelastic  Fluid  Law  For  Ice 
L.  W.  Morland  and  U.  Spring 

School  of  Mathematics  and  Physics 
University  of  East  Anglia,  Norwich,  (U.K.) 

Abstract 

A  frame-indifferent  differential  operator  law  relating 
stress,  stress-rate,  strain-rate,  and  strain-acceleration,  is 
constructed  to  describe  the  qualitative  features  of  both  constant 
stress  and  constant  strain-rate  response  in  uni-axial  stress 
experiments.  The  structure  of  the  tensor  relation  cannot  be 
determined  by  uni-axial  data,  and  a  variety  of  models  reducing 
to  the  required  uni-axial  form  are  presented.  The  differential 
law  allows  exact  description  of  a  family  of  strain-rate  curves 
at  different  constant  uni-axial  stress,  and  some  features  of 
the  stress  curves  at  different  constant  strain-rates,  showing 
that  constant  stress  and  constant  strain-rate  response  reflect 
some  common  material  properties,  but  also  reflect  independent 
properties.  That  is,  both  types  of  response  are  necessary  to 
define  the  material.  A  reduced  model  which  adopts  a  restricted 
form  of  dependence  on  stress  and  strain-rate  invariants  is  also 
analysed.  Idealised  families  of  constant-stress  responses 
are  constructed  and  various  forms  of  limited  matching  to  the 
reduced  model  are  presented  to  show  the  extent  to  which  the 
simplified  law  can  describe  the  complete  response. 
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Introduction 

In  ice  engineering  the  mechanical  response  of  ice  under 

maintained  stress  is  required  over  times  sufficient  for  creep 

to  strains  of  a  few  percent,  or  to  rupture  at  appropriate 

stress  combinations.  This  covers  the  instantaneous  elastic 

strain,  primary  decelerating  creep,  secondary  or  approximately 

steady  creep,  and  some  part  of  the  accelerating  tertiary 

creep  which  on  glaciological  time  scales  is  supposed  to 

approach  steady  creep.  The  Young's  modulus  and  bulk  modulus 

lO  -  7 

of  ice  are  of  order  10  Nm  (Sinha  1978,  Mellor  1980),  so 

6  -2 

that  a  moderate  uni-axial  compressive  stress  of  10  Nm 

-4 

(10  bars)  induces  an  elastic  strain  of  only  10  compared 

-2 

with  creep  strains  of  10  common  in  application,  and  it 

is  convenient  to  neglect  the  elastic  response  and  model  only 

the  creep.  The  constitutive  law  may  still  be  of  "solid  type", 

depending  on  a  reference  configuration,  which  is  essential 

to  model  anisotropy.  While  both  initial  and  induced  anistropy 

can  be  significant  there  has  been  no  systematic  investigation 

(Mellor  1980) ,  so  a  realistic  constitutive  model  is  not  yet 

feasible.  Here  we  focus  on  a  differential  operator  "fluid  type" 

model  to  describe  observed  viscoelastic  creep,  necessarily 

isotropic  in  all  configurations. 

Constant  uni-axial  stress  experiments  to  determine  the 

short  time,  small  strain,  creep  have  been  reported  by  Sinha 

(1978a,  1978b)  and  Gold  and  Sinha  (1980) .  Total  creep  times 

were  typically  of  several  minutes  duration  with  maximum 

-4 

strains  of  order  3  *  10  , 


but  a  wide  range  of  constant 
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temperature,  228K  to  263K,  was  covered.  The  results  are 
used  to  infer  an  empirical  relation  for  the  uni-axial  strain 
at  time  t  in  terms  of  the  stress  and  time  explicitly,  which 
incorporates  instantaneous  elastic  response,  recoverable  creep, 
and  viscous  flow.  It  is  not  shown  how  the  relation  can  be 
related  to  a  three-dimensional  tensor  law  (coordinate  invariant) 
necessary  for  combined  stress  loading,  nor  is  it  clear  how 
the  response  to  a  non-constant  stress  history  is  deduced  since 
there  is  no  superposition  principle  for  non-linear  response. 

We  will  now  examine  a  frame  indifferent  tensor  relation  which 
contains  the  minimum  dependence  on  rates  necessary  to  describe 
the  main  qualitative  features  of  uni-axial  response.  As 
demonstrated  in  an  earlier  more  limited  treatment  (Morland  1979), 
uni-axial  response  cannot  determine  the  shape  of  a  tensor 
relation,  which  would  require  tri-axial  (three  independent 
principal  stresses),  or  combined  shear  and  axial  stress,  tests. 
There  are  an  infinity  of  tensor  relations  which  reduce  to  the 
same  uni-axial  response,  so  that  construction  of  an  explicit 
tensor  relation  requires  many  more  assumptions  about  shape  and 
the  arguments  of  response  functions.  The  various  restrictions 
and  possible  alternatives  are  noted  in  our  illustrative 
construction. 

Mellor  (1980)  has  described  qualitatively  the  uni-axial 
response  at  constant  temperature  under  constant  compressive 
stress  and  under  constant  compression  strain-rate.  He  notes 
that  data  from  both  tests  must  reflect  the  same  material 
properties,  but  points  out  that  no  common  rheological  model 
simulates  the  observed  responses  in  both  tests.  We  now  construct 
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a  differential  operator  viscoelastic  fluid  law  which  describes 
the  main  features  of  the  responses  in  both  tests,  showing 
that  neither  the  constant  stress  nor  constant  strain-rate 
data  can  determine  the  law  completely,  so  that  these  two  forms 
of  response  do  indeed  contain  independent  information  and 
neither  reflects  fully  the  mechanical  properties.  This  contrasts 
with  the  theory  of  linear  viscoelasticity  for  which  the 
existence  of  a  linear  hereditary  integral  law  (Gurtin  and 
Sternberg,  1962)  implies  that  a  unique  description  is  determined 
by  the  creep  functions,  namely  response  under  constant  stress. 

Figure  1  shows  the  axial  strain-rate  r(t)  for  a 

uni-axial  compressive  stress  o  applied  at  time  t  =  0  and 

maintained  constant.  In  the  primary  creep  am  ,  r  decreases 

from  the  initial  strain-rate  rQ(o)  to  the  minimum  strain-rate 

rm(a),  both,  in  general,  depending  on  the  stress  level  a  , 

and  the  turning  point  m  corresponds  to  the  inflexion  point 

(or  secondary  creep)  on  a  strain-time  curve.  In  the  tertiary 

creep  me,  r  increases  from  rm(o),  and  we  suppose  an 

equilibrium  strain-rate  re(a)  is  attained  asymptotically, 

corresponding  to  the  commonly  assumed  steady  creep  appropriate 

to  many  glaciology  applications.  In  Fig.  1  we  have  shown 

re(o)  <  rQ{0),  and  will  present  analysis  for  this  case,  but 

a  similar  treatment  can  be  given  for  the  case  re(a)  >  rD(o). 

The  time  to  minimum  strain-rate  is  t  (o),  which  increases 

m 

significantly  as  a  -*■  0  so  that  minimum  strain-rate  is  not 
attained  in  many  low  stress  laboratory  tests.  A  steady  creep 
law  for  glaciology  applications  should  be  the  tertiary  limit 


* 
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r  (o),  not  the  minimum  strain-rate  r  (a) ,  but  Mellor  (1980) 
e  m 

suggests  that  the  distinction  diminishes  as  stress  decreases 

and  may  not  be  significant  at  deviatoric  stresses  below 
5  -2 

10  Nm  (1  bar)  typical  of  many  glaciology  situations. 

Figure  2  shows  the  uni-axial  compressive  stress  a(t) 
when  the  axial  strain-rate  r  is  held  constant.  In  the  first 
stage  OM  ,  o  rises  from  zero  to  a  maximum  stress  oM(r) 
at  time  tM(r)  •  then  decreases  in  the  stage  ME,  to  an 
asymptotic  equilibrium  limit  0  (r)  for  consistency  with 
a  tertiary  creep  limit  rg(0).  Mellor  (1980)  also  suggests 
that  aE  ■+  oM  as  strain-rate  •*>  0.  The  assumption  of  a 
steady  (finite  rate)  creep  re(o)  as  t  -*•  00  at  constant 
stress  0  implies  a  correspondence  with  the  limit  o£(r) 
at  constant  strain-rate  r  given  by 

re[°E(r)J  “  r  and  oErre(a)l  =  a  , 

or  r”1^)  =  0£(  r)  and  ©“^(o)  =  re(o). 

Mellor  (1980)  suggests  that  there  are  indications  that  the 
maximum  stress  at  a  given  constant  strain-rate  is  the  constant 
stress  required  to  produce  that  strain-rate  as  the  minimum. 

We  will  adopt  this  property,  thus 

rm&Mlr)]  *  r  and  oM|>m<a)]  -  o  .  (j) 

°r  r„  (r>  =  dM(r)  an<3  aM1(0)  '  rm(°)  * 

m  M  Mm 

Both  inverse  relations  (1)  and  (2)  are  illustrated  in  Fig.  3. 

Mellor  further  speculates  that  the  strain  e  (0)  at  the 

m 

minimum  strain-rate  in  the  constant  stress  test,  and  the  strain 
c^(r)  at  the  maximum  stress  in  the  constant  strain-rate  test. 


do  not  vary  significantly  with  a  and  r  respectively  over 
some  range  of  a  and  r  ,  and  are  approximately  the  same, 
slightly  less  than  0*01. 

While  corresponding  tensile  tests  are  not  practical 
over  a  wide  stress  range  it  is  expected  that  tensile  response 
will  differ  in  magnitude,  and  be  limited  by  tensile  fracture. 

In  our  simple  illustrative  model  we  restrict  correlation 
to  the  compression  responses  shown  by  Figs  1  and  2,  but 
indicate  how  distinct  independent  compression  and  tension 
response  (if  known)  will  influence  the  construction. 

A  viscous  fluid  model  prescribes  a  as  a  function  of 
r  ,  or  vice-versa,  and  implies  constant  r  for  constant  a  , 
and  constant  o  for  constant  r  ,  in  contrast  to  the  variation 
with  time  shown  in  Figs  1  and  2.  By  allowing  stress  to 
depend  also  on  strain-acceleration ,  so  that  a  is  given  as 
a  function  of  r  and  r  ,  the  primary  creep  am  in  Fig.  1 
of  a  constant  stress  test  can  be  modelled,  as  demonstrated 
by  Morland  (1979)  in  the  case  when  the  minimum  strain-rate  is 
an  asymptotic  limit;  that  is,  accelerating  tertiary  creep 
is  absent.  However,  we  will  show  that  such  a  first  order 
differential  relation  for  r  in  terms  of  o  can  describe 
the  full  constant  stress  response  of  Fig.  1.  In  contrast, 
the  stress  variation  at  constant  strain-rate  shown  in  Fig.  2 
cannot  be  described  by  such  a  model  which  reduces  simply  to 
a  viscous  fluid  law,  o  a  function  of  r  ,  when  r  -  0.  By 
analogy  with  the  transient  creep  construction,  it  is  necessary 
to  incorporate  dependence  also  on  stress-rate,  so  that  when 
r  =  0  the  reduced  form  is  a  first  order  differential  relation 


for  o  in  terms  of  r  .  We  will  now  construct  a  frame 


indifferent  tensor  relation  relating  stress,  stress-rate, 
strain-rate,  strain-acceleration  which  can  model  both  families 
of  response  curves,  and  which  necessarily  has  the  minimal 
rate-dependence  for  this  purpose.  The  complete  constant  stress 
response  can  be  matched,  but  completion  of  the  model 

requires  some,  but  not  all,  of  the  constant  strain-rate 
response,  or  vice-versa.  That  is,  constant  stress  and  constant 
strain-rate  responses  reflect  some,  but  not  all,  of  the  material 
properties.  We  also  consider  a  simpler  reduced  form,  and 
demonstrate  by  illustrations  the  extent  to  which  it  reproduces 
a  set  of  idealised  responses. 

Differential  operator  law 

The  conventional  incompressibility  approximation  is 
made  so  that  the  Cauchy  stress  a  is  determined  by  the 
deformation  history  only  within  an  arbitrary  additive  isotropic 
pressure  p  1  .  That  is,  the  constitutive  law  determines  only 
the  deviatoric  stress 

S  =  a  -  j  tra  1  ,  tr  S  «  0  .  (3) 

The  general  viscous  fluid  law  can  be  written  (Morland  1979) 

5  =  MV?  +  MV^?2  “  §I2il  '  <4> 

or  equivalently 

?  =  MV?  +  MV^?2  ■  §J2$  *  (5) 

where  the  strain  rate  tensor  D  is  the  symmetric  part  of  the 
spatial  velocity  gradient,  and  the  non-trivial  strain-rate 
invariants  1^  and  deviatoric  stress  invariants  (k  =  2,3) 


^ _ 
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are 

I2  =  ■!  tr  D2  ,  I3  =  det  D  ,  J2  =  |  tr  S2 ,  J3  =  det  S  .  (6) 

The  forms  (4),  (5)  are  consistent  with  the  requirements 
=  tr  S  =  0  ,  1^  =  tr  D  =  0  (incompressibility) ,  and  the 

response  coefficients  <J>2  and  i|^,  \|>2  depend  only  on 

two  invariants  I2,  I3,  and  J2,  respectively.  Analytic 

inversion  of  the  general  forms  (4),  (5)  is  not  possible,  but  a 
conventional  glaciology  flow  law  assumes  ij>2  =0,  ij;^  =  <|>^(J2), 
which  implies  <J>2  =  0,  4> ^  ^  ( X2 )  given  implicitly  by 

<t>i <i2)  =  1-  (?) 

Clearly,  uni-axial  response  cannot  determine  two  functions 
<J>i#  <f>2,  or  ij^,  ij»2,  each  depending  on  two  arguments  I2,  I3, 
or  J2#  J3,  respectively , and  furthermore,  bi-axial  tests 
determine  only  the  same  combination  (3I2)^$^  +  I2<{>2  on 
I3  =  ±2(I2/3)^  ,  or  (3J2)^i|»^  +  J2^2  °n  J3  =  ±2(J2/3)7  ,  with 
two  independent  principal  stresses  a^,  o2  =  (Morland  1979). 

This  is  a  direct  consequence  of  incompressibility  tr  D  =  0  which 
imposes  one  relation  between  the  principal  strain-rates, 
d^  +  2d2  =0  in  bi-axial  deformation,  and  correspondingly 
S1  +  ^s2  =  ®  from  tr  S  =  0  ,  so  that  each  principal  component  of 
the  law  (4)  or  (5)  is  a  multiple  of  the  others.  Thus,  tri-axial 
tests  with  three  independent  principal  stress  components  a2, 

o3,  or  combined  axial  and  shear  stress  tests  (shear  stress  and 
shear  strain-rate  do  not  enter  tr  S  and  tr  D ) ,  are  required  to 
determine  two  functions  of  two  arguments.  This  lack  of  determinancy 


of  tensor  relation  shape  by  uni-axial  data  is  compounded 
when  strain-acceleration  and  stress  rates  are  incorporated 
in  the  law. 

A  varying  strain-rate  r(t)  under  constant  stress, 

Fig.  1,  can  be  modelled  by  including  stress  dependence  on 
the  strain-acceleration  of  a  material  element,  measured  by 
the  second  order  Rivlin-Ericksen  tensor 

A(2^  =  2D  +  4D2  +  2 (D  W  -  W  D)  ,  (8) 

where  the  rotation  rate  W  is  the  skew  part  of  the  spatial 
velocity  gradient 

L  =  D  +  W  .  (9) 

(2) 

A  is  a  frame  indifferent  tensor.  Note  that 

tr  A(2)  =  4 tr  D2  *  8l2  (10) 

since  tr D  =  0,  tr(D  W)  =  tr(W  D) ,  while  second  and  third 

invariants  will  depend  on  i2,  i-j.  All  products  of  powers  of 

( 2 ) 

D  with  powers  of  A  ,  with  coefficients  depending  on 

(2 ) 

invariants  of  D,  A  and  such  products,  are  frame  indifferent, 

so  there  is  no  simple  compact  frame  indifferent  expression. 

Following  Morland  (1979),  the  dependence  on  strain-acceleration 

(2) 

will  be  limited  to  the  first  power  A  ,  but  here,  the 
coefficient  is  allowed  to  depend  on  i2  as  well  as  on  I2,  I^, 

J2»  J3/  to  construct  one  model  which  describes  both  primary 
and  tertiary  creep. 

Similarly,  a  varying  stress  o(t)  at  constant  strain-rate, 
Fig.  2,  can  be  modelled  by  including  a  frame  indifferent  deviatoric 
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stress  rate 

S(1)  =  S  +  S(D  +  W)  +  (D  -  W)S,  trS(1)  =  2tr(S  D)  ,  (11 

with  coefficient  a  function  of  all  invariants,  but  here  limited 
to  I 2,  I3,  J 2,  J3,  J2.  Thus,  an  appropriate  frame  indifferent 
relation  between  D,  A^,  S,  S^,  is 

+  ^3[s(1)  -  §tr(s  D)l]  =  +  4>2[D2  -  §I2J;]  +  4>3[D  +  DW-  W  D]  (12 

2  ( 2 ) 

where  the  D  term  of  A  is  incorporated  in  the  <f>2  term.  Note 
that  A(2)  =  0  *  6  =  0,  even  for  W  =  0.  In  one  model  4^,  <J>2, 
t|/^  are  functions  of  I2,  I3,  J2,  J3,  but  <J>3,  ip^  depend  also 
on  i2,  J2,  respectively.  An  alternative  dependence  on 
strain-acceleration  and  stress-rate  only  through  the  rate  invariants 
i2,  J2,  is  obtained  by  setting  <j>3  =  if>3  =  0  and  allowing 
\p^  to  depend  on  the  rate  invariants.  The  mixed  combinations 
ip  j  =  0,  <(>3  ?  0,  and  $3  =  0,  \p^  /  0  are  also  considered. 

Again,  uni-axial  response  cannot  determine  such  a  shape 
distinction. 

Uni-axial  response 

Now  consider  a  uni-axial  compressive  stress  -o^  =  o  >  0, 
other  cKj  =  0,  with  corresponding  strain-rate  -D-^  =  r  >  0, 

°22  =  °33  =  '  other  Dij  =  and  zero  rotation-rate  W  =  0. 

Then 


•  -* 

1 

O  N> 
Q 

0 

0 

0 

0 

,  s(1)  =  s  +  i 

4ar 

0 

0 

or 

0 

0 

i  0 

0 

a 

0 

0 

or 
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? =  1 


-2r 

O 

O 

2  1 

4r  0  0  > 

2 

0 

r  0 

'  ?  =  T 

0  r  0 

2 

,  0 

0  r 

1 

u 

0 

0 

(2)  _ 


=  2D  +  4D  , (14 


-r  -  1  2  _  _  2  3  _  _  -  ,1  7  >2  •  _  2  •  c(1>  ~ 

J 2  ~  o  i  J3  -  -  27  0  2(3^2)  /  J2  ^  o  (j  ;  tr  S  •*  2r  o  1  1*5 


3  2' 
3 


_  _  3  2  T_  1  3  _  _  ,1_  .2  •  _  3  • 

*2  —  4^  *  ^3  ~  —  4r  ~  ”  2  3^2  '  7  2  2r  r  * 


(16 


Each  principal  component  of  the  law  (12)  gives 

f  ♦l0  +  §  ^3(0  ~  ra  )  =  ^ir  “  \  <t>2r2  +  <f>3*‘  d7 

where  <J> ^ ,  <J>2,  <J>3  depend  on  r  and  r  through  (16),  and  ,  <(/3 


depend  on  a,  a, 

and  r 

through  (15). 

For  constant 

stress. 

o=0,  (17)  reduces  to 

4>3r  +  ^r 

-  W2 

+  §*3  o  r  =  o  , 

(18 

which  is  a  first  order  differential  equation  for  r(t)  given  the 
constant  a  and  an  initial  condition 

r (0+)  =  rQ(o)  .  (19 

Similarly  for  constant  strain-rate,  r  =  0,  (17)  reduces  to 

§  <1» 3O  +  •f'J'jO  -  §^3ro  *  $ir  "  ^2r2  (2° 

which  is  a  first  order  differential  equation  for  o(t)  given  the 


constant  r  and  an  initial  condition 


a (0+)  =  0.  (21) 

The  non-linear  differential  equations  (18)  and  (20)  must  have 
unique  positive  solutions  r(t),  a(t)  in  t  >  0,  and  for 
each  constant  o  and  constant  r  respectively,  describe  qualit¬ 
atively  at  least  a  class  of  response  curves  of  the  form  shown  in 
Fig.  1  and  a  class  of  response  curves  of  the  form  shown  in  Fig.  2. 
While  it  was  not  anticipated  that  a  first  order  differential  law 
would  match  a  given  family  of  constant  stress  responses  exactly,  we 
find  this  can  be  achieved  by  appropriate  choice  of  the  coefficients 
in  (18)  .  In  addition,  (20)  can  then  match  the  family  of  primary 

stress  increases  at  constant  strain-rates,  but  predicts  the  stress 
relaxation.  ,  .  . 

A  symmetric  response  in  tension  requires  o  ,  a  to  change 
sign  with  r,  r  ,  which  is  a  restriction  on  the  a,  r  dependence 
of  the  response  functions.  Independent  compression  and  tension 
response  will  in  general  require  dependence  on  both  even  and 
odd  invariant^  in  the  uni-axial  reductions  (15),  (16),  where 

J3  =  ^  an<*  *3  =  +2^1I2  ^  for  tension-  $1'  ♦ 2 '  ^3' 

determined  as  functions  of  r,  a  by  correlation  with 

compression  data  can,  in  principle,  be  expressed  as  functions  of 

I2  or  Ij»  or  a  combination,  and  J2  or  Jj,  or  a  combination, 

but  different  dependences  on  the  invariants  will  predict  different 

tension  responses.  Hence,  complete  knowledge  of  both  tension 

and  compression  response  provides  a  restriction  on  I2,  I3  and 

J 2#  J-j  dependence.  A  simple  example  is  the  viscous  law  (4), 

2  12 

which  becomes  3-0  =  -  ■j$2r  •  subject  to  symmetric  tension- 

compression  response  requiring  o/r  invariant.  Then  the 


J 


I. 


combination  -  ^r4>2  must  be  invariant  as  r  changes  sign, 

2 

suggesting  that  <p^  depends  on  I2  and  I3  and  <j>2  is  an 

2 

odd  power  of  I3  multiplied  by  a  function  of  I2  and  I3. 


This  assumes  smooth  response  as  r  -*■  ±0.  We  now  focus  on  the 
laws  (18)  and  (20)for  uni-axial  compression,  o  >  0,  r  >  0. 


Constant  stress  response 


An  essential  feature  of  the  creep  curve  shown  in  Fig.  1 

is  the  repeat  of  the  strain-rate  values  in  r  (o)  <  r  <  r  (o) 

me 

during  tertiary  creep  t  >  tm(o).  If  we  considered  re(°)  >  rQ(a), 
then  the  duplicated  values  would  be  in  the  primary  creep  range 
r  (o)  >  r  >  r  (cr)  .  Thus,  the  differential  equation  (18)  must 
describe  a  double-valued  r  for  rm(°)  <  r  <  re(o),  with 
r  =  r_  <  0  during  the  primary  creep,  t  <  tm ( a ) ,  and 
r  =  r+  >  0  during  the  tertiary  creep  t  >  tm(o),  and  r ( tm)  -  0. 
The  necessary  form  compatible  with  (18)  is 


r2  +  f (r,o) r  =  F(r,a) , 
which  implies 

2r  =  -f  i  {f2  +  4F)^  , 


r  (0)  =  rQ(a)  >  0, 


(22) 


(23) 


The  relation  of  f  and  F  to  the  response  coefficients  is 
deferred  until  the  constant  strain-rate  equation  has  also  been 
treated. 


By  (23)  we  see  that  for  each  o  and  each  r  that  F  >  0 
is  a  necessary  and  sufficient  condition  for  one  root  t_  <  0  and 
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one  root  r+  >  0;  hence 

F(r ,o)  >  0,  rm (a )  <  r  <  re(o)  .  (24) 

For  r__  -*■  0  and  r+  -*■  0  as  r  ■*  r^fa)  ,  it  is  necessary  that 

frrm(o),  cr]  =  F[rm(o),o]  =  0  ,  (25) 

and  r  switches  from  r  to  r.  as  t  passes  t  (a) 

-  +  m 

provided  that 

8r  =  [l  +  f(f2  +  4  F )  ”  ^  J  1 2  ||(f2  +  4F)*5  +  ^(f2  +  4F)j 

>0  as  t  -*■  t  (o)-  ;  (26) 

m 

that  is,  if 

3f/3r  >  O,  3F/3r  >0,  at  r  =  r  (o)  .  (27) 

Since  f(r,o)  and  F(r,o)  have  no  role  in  r  <  r  (o),  it 

m 

is  convenient  to  extend  them  continuously  by 

f (r,o)  =  F (r , o)  =  O,  r  «  rm(o).  (28) 

The  asymptotic  condition  r  rg(o)  as  t  -*•  »  requires 
r+  •*  O  as  r  -*•  re(o),  which  implies 

F[r#(0),0]  =  0,  (29) 

and  since  there  is  only  primary  response  in  r  >  r Q(a)  we  set 

F(r ,0)  HO,  r  >  rg(o).  (30) 

By  definition,  the  time  to  approach  re(o)  in  tertiary  creep  is 


unbounded;  that  is 


V01 

r 

2dr ' 

as  r  +  r 

e 

i 

:m(0> 

[f2(r',0) 

+ 

4F  (r  * ,  0  )J  ** 

-  f  (r  ’ 

,0) 

while 

the  time 

V°> 

to 

minimum  strain-rate  is 

bounded ; 

ro<0) 

t(o) 

in 

[ 

2dr ' 

_  < 

i 

1 

r_(0) 

xn 

[f2(r  *,0) 

+ 

4F(r '  ,0)]  ** 

+  f  (r ' 

,0) 

00  • 

(3 


Given  that  fCre(a)»a3  =  -  r_  at  r  =  r e(o)  is  not  zero  and 
f  is  analytic  there,  the  unbounded  integral  (31)  requires 


F(r,cr)  „<  0(r  -  r)  as  r  -»•  r  (cr), 

6  e 


(3; 


and  the  bounded  integral  (32)  requires 

Hr, a)  +  [f2(r ,o)  +  4F(r,a)J!5  >  0( r  -  rm>  as  r  >  rm(0),  (3- 

together  satisfied  by 


f  *  0(r  -  r  ) ,  F  =  0(r  -  r  )  as  r  r 

m  m  m, 

F  =  0(r  -  r)  as  r  -*■  r  . 

e  e 

Note  that  condition  (30)  implies  the  reduction  of  the 
differential  equation  (22)  on  rQ(0)  >  r  >  re(o)  to 

r  Qr  +  f(r,o)3  =  0,  r(0)  =  rQ(o), 


which  permits  a  solution  r  =  ro(o).  The  branch  =  -  f(r,o), 
with  f  >  0,  determines  the  primary  creep.  At  the  branch  switch 
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tm(o>,  F  *  O,  r+  =  r_  =  O,  but  r  >  0,  and  r+  becomes 
positive. 

For  numerical  correlation  with  an  idealised  family  of 
uni-axial  response  we  consider  a  reduced  model  involving 
functions  of  single  argument,  which  could  be  expressed  in  terms 
of  invariants  I2,  J2-  Define 

<*1  -  r2  -  (*2  *  1>2  -  rm(o0  l>e<°>  ‘  r2]  •  <37> 

and  restrict  r  dependence  by 

f(r,o)  =  q  (o )  1(0^)  =  0  for  «  0, 

(38) 

F(r,o)  =  q2 (o )  F  (Qg),  F(^2)  =  0  for  Q2  «  O, 

which  satisfy  the  requirements  (25),  (29),  and  subsidiary 
conditions  (28),  (30),  identically.  The  choice  of  arguments 
Q2'  is  for  direct  expression  in  terms  of  l2>  The 
conditions  (35)  become 

F(02)  =  0(fl2)  as  Q2  -*■  0,  f(G1)  =  0(0^  as  -*•  0.  (39) 

As  r  increases  from  r^  to  r*  =  [^(r2  +  r2)]^,  0l2  increases 

12  2 

from  zero  to  ^(re  -  rm) ,  then  decreases  to  zero  as  r  increases 

from  r*  to  rg.  Thus  ^2(r)  has  duplicated  values  in 

r  4  r  *  r  ,  which  will  restrict  correlation  over  the  full 
m  e 

primary- tertiary  creep  response.  As  r  increases  from  r  to 

in 

rQ,  is  strictly  increasing.  With  the  forms  (38), 


1 
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ro(o) 


t_(o) 

ID 


2dr ' 


1  f  _ _ 

^  Jr  (0)  (f2(ai)  +  4F(a^)]!5  +  f(d{] ' 


(40) 


Gi  -  r'2  -  rm  '  Gj  -  (r'2  -  r2)  (r2  -  r-h 


Constant  strain-rate  response 

Analogous  to  the  constant  stress  analysis,  the  necessary 
form  compatible  with  (20)  to  describe  o (t)  at  constant  r 
shown  in  Fig.  2  is 

*  2  $  • 
a  -  g (r,o)o  =  G(r , 

which  implies 

2a  -  g  ±  (g2  +  4G) * , 

For  a+  >  0  and  o_  <  0 
require 

G(r,a)  >  0,  aE(r)  <  °  *  (43) 

while  a+  ■*  0  and  o_  0  as  a  -*■  aM(r)  requires 

g[r,aM(r)']  =  G[r,oM(r)]  =  0.  (44) 

a  switches  from  o+  to  a_  as  t  passes  through  tM(r) 
provided  that  a  <  0  there*,  that  is,  if 

If  >  0  ,  If  >  0  at  o  -  0M(r>.  (45) 


a) ,  a  (0)  =  0  , 


(41) 


-  2(  M  • 


(42) 


on  the  duplicated  stress  range  we 


The  asymptotic  condition  a  oE(r)  as  t  -*■  00  requires 
o_  -*•  0  as  a  oE(r),  which  implies 

G[r,oE(r)]  =  0  , 

and  an  infinite  relaxation  time  requires,  given 
g|lr,oE(r)J  =  o+  at  c  =  oE(r)  is  not  zero  and  g  is 
analytic  there, 


G(r,o)  ^  0  (o  -  oE)  as  a  -*•  o£(r). 


The  time  to  maximum  stress  is 
o„(r) 


r)  m  f  M  _ 2da 1  _ 

^  r  Jq  [g 2(r^a')  +  4G(r,a')]*  +  g(r,o') 


provided  that 


g(r,o)  +  jjg2(r,a)  +  4G(r,a)]^  >  0{oM  -  o)  as  a  oM(r) 


Conditions  (47)  and  (49)  are  satisfied  by 


g  =  0(oM  -  o),  G  ■  0(cM  ”  o)  as  o  +  oM, 


G  *  0(o  -  o£)  as  o  -*■  o£ 


We  can  set 


g(r,o)  =  G (r,o)  =  O,  o  *  oM(r), 


G(r,o)  =0,  a  4  o£(r) . 


Mow  on  0  404  o£(r),  (41)  becomes 

o[o  -  g(r,o)J  *  0,  o(0)  »  O, 
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which  allows  a  solution  a  =  0.  However,  the  primary  response 

is  qiven  by  a.  =  g(r,o),  and  at  the  branch  switch  t,,(r) 

•  •  •* 
where  G  =  0,  o  =  0,  a  becomes  negative  since  o  <  0. 

Now  we  apply  the  inverse  relations  (1)  and  (2)  which  state 

that  the  curves  r  =  r  (o)  and  o  =  a„(r)  in  an  (r,o)  plane 

are  identical,  and  the  curves  r  =  r  (a)  and  a  =  o^tr)  are 

G  E 

identical.  Furthermore,  r  <  rm(a),  or  <  0,  corresponds 

to  a  >  aM(r),  and  r  >  re(a),  or  CL^  >  O  and  <  O, 
corresponds  to  o  <  oE(r).  Thus  the  conditions  (51)  become 


g(r,a)  5  G(r,a)  50,  r  <  rm(a), 
G(r,o)  5  0,  r  >,  re(o)  . 


(53) 


Conditions  (28),  (30),  (53)  are  illustrated  in  Fig.  3. 

Corresponding  to  the  reduced  model  (36),  (37)  we  define 


g  (r,  a)  =  s  (r)q(o)g  (C^)  ,  g^)  5  0  for  0, 

G(r,a)  =  q2  (a ) G (#2 )  ,  G (Q^)  =  0  for  ^  <c  0,  (54) 

G(fl2)  =  0(a2)  as  Q2  -  0,  gf^)  =  O^)  as  -  O, 

which  satisfy  the  requirements  (53)  identically,  and  provided 
that  r^fcr),  o^(r)  $  0,  satisfy  condition  (50).  It  is  found 
that  the  extra  factor  s(r),  here,  or  as  an  f  factor  in  (38), 
is  necessary  for  time  scale  control. 


Response  coefficients 

We  now  construct  four  sets  of  response  coefficients  ^ ^ , 

♦l#  $2'  ^3  which  reduce  the  constant  stress  and  constant  strain-rate 
relations  (18)  and  (20)  to  the  required  forms  (22)  and  (41) .  The 
first  supposes  that  ^  /  0  so  that  there  is  rate  dependence 
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through  both  tensors  and  S^,  and  the  second 

assumes  <J>3  =  4/^  =  0  so  that  rate  dependence  is  through 
only.  The  other  two  are  the  combinations  4*3  =0,  $3  ¥  0, 
and  =  0,  ^3  ^  0.  With  each  set  there  is  the  implication 

F(r,a)  =  G(r,a),  shown  in  Fig.  3,  which  represents  the  coupling 
of  constant  stress  and  constant  strain-rate  responses  in  this 
model,  or  the  extent  to  which  they  hinge  on  a  common  material 
property. 

In  the  case  4> ?  0,  4*3  7*  0  we  can  take  4/^  =  1  without 

loss  of  generality.  Then,  restricting  rate  dependence  to  the 
new  <p 2 >  <P 3  and  comparing  (18)  and  (22),  shows  that 

<P3  =  r  +  f(r,o),  =  f(r,o)  when  r  =  0,  (55) 

so  <f>3  is  a  function  of  1^,  I2,  J2 ,  I3,  J3  in  general. 
Comparing  (20)  and  (40)  we  find  4*3  is  a  function  of 
J2,  I2,  J2,  I3,  J3  given  by 

2 

-j4>3  =  -  o  +  g(r,a)  -  r  o  ,  =  g(r,a)  -  ra  when  0  =  0  ;  (56) 

this  sign  choice  gives  G  =  +  F.  Then  (2) ,  (41)  are  recovered 
with 

F  (r ,  o )  =  G  (r ,  0 )  =  -  tj>  x  r  +  i(f>2r2  “  rag(r,a)  +  r2o2  +  ~  a  ,  (57) 

so  that  <l>^,  4> 2  are  functions  of  I2,  J2,  I^,  only.  Again, 

only  a  combination  of  <J>^,  4>2  is  given  by  F. 

In  the  case  <J>3  =  ^3  =  0,  restricting  rate  dependence  to 
4>^,  the  above  comparisons  show  that 

0  =  -  r2  -  f  (r,a)  r  -  a2  +  g(r  ,a)o  +  CUr,a), 

J  1  (58) 

F(r,a)  =  G(r,o)  =  -  ♦j^r  +  ^2r2  +  Ct(r,a); 
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where  depends  on  i2,  J2,  I2'  J2'  I3'  J3'  and  a  com^ination 

of  4*^  and  <f>2  depends  on  I2,  J2,  l3,  J3>  Normalising  by 

4^  =  1  changes  <j>^,  4>2  but  reproduces  the  tensor  law  (12) 

with  coefficients  given  by  (58).  The  arbitrary  0L(r,a) 

represents  in  (12)  an  S  term  with  rate-independent  coefficient, 

2 

with  corresponding  D,  D  terms.  This  generality  is  lost  in 
the  first  construction  where  rate-dependence  is  restricted  to 
<J>3,  4>3  after  setting  4^  =  1.  Clearly,  the  uni-axial  response 
functions  f,  g,  F  cannot  determine  #,(r ,o)  in  this  model. 
Alternative  models  with  4>2  rate  dependent  can  be  constructed 

Similarly,  again  excluding  rate-dependence  from  4^,  <j>2> 
for  <J>3  5^  O,  ip3  =  0, 

4>3  =  r  +  f  (r,o)  , 

=  -  o2  +  g  (xfa)b  +  b(r,o)  ,  (59) 

F(r,o)  =  G(r,o )  =  -  <f>1r  +  j<P2r2  +  t>(r,c), 

and  for  <J>3  =  0,  /  O, 

2  . 

3  <ji3  =  -  a  +  g  (r,o )  -  r  a  , 

=  -  r2  -  f(r,a)r  +  c(r,a),  (60) 

F (r,o)  =  G(r,a)  =  -  (f^r  +  i<fr2r2  +  r2cr2  "  rag(r,a)  +  c(r,a) 

where  b(r,a),  c(r,a)  are  arbitrary. 

The  variety  of  dependence  of  the  response  coefficients  on 
r  and  a  eliminates  simple  dependence  on  integral  powers  of  I2, 
J2,  but  including  I3,  J3  allows  rational  functions  of  integral 
powers.  Dependence  on  tr(S  D)  =  ro  could  also  be  incorporated. 
Uni-axial  response  does  not  distinguish  the  various  forms  of 
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dependence.  Each  of  the  four  models  described  above  will  predict 
different  response  in  any  non-uni-axial  stress  geometry,  for 
example,  simple  shear  stress  or  simple  shear  motion,  and  different 
dependence  on  invariants  within  each  model  similarly  predicts 
different  response.  Finally,  following  Morland  (1979),  a 
tentative  initial  condition  for  the  tensor  law  (12)  when  stress 
is  prescribed  is 

.  ,  ,  r  (o) 

D(0)  =  h(J2)S(0),  h(±c  )  =  j  -  (61) 

Uni-axial  response  functions 

It  is  evident  that  three  functions  f(r,o),  g(r,o),  F(r,a) 
cannot  exactly  reproduce  independent  strain-rate  curves  r(t) 
for  each  constant  a  and  stress  curves  a(t)  for  each  constant 
r.  A  differential  operator  law  including  only  strain-acceleration 
and  strain-rate  evaluated  at  current  time  is  a  restricted  history 
dependence.  We  can  incorporate  the  important  properties  rm(a), 
rg(a)  through  the  Q.^,  definitions  (37),  and  the  associated 

properties  oM(r),  oE(a)  through  the  inverse  relations  (1). 
Because  of  the  restriction  F  =  G  there  are  two  functions  f,  F, 
available  to  correlate  with  constant  stress  response,  but  then  only 
one  function  is  available  to  correlate  with  constant  strain-rate 
response,  or  vice-versa.  Thus  one  class  of  response  is  given 
priority.  We  will  develop  the  former  case,  but  a  similar 
procedure  applies  to  the  latter. 

For  each  constant  o  the  given  strain-rate  response 
determines  r+(t)  =  R+(r,o)  >,  0  and  r_(t)  =  R_(r,o)  0.  By 

(23), 

rm  ■»<  r  4  rQ:  f  *  -(R+  +  R_),  (f2  +  4F)*  =  R+  -  R 


r 


(62) 
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with  F  =  0  and  R+  =  O  for  re  ^  r  ^  rQ.  Recalling  the 
sufficient  analytic  requirement  (35),  we  can  write 

f  =  0]f  *  F  =  CL2F  .  (63) 

A  A 

where  f,  F  are  non-zero  at  r  =  rm(o)  and  r  =  re(c),  and 

factor  &,  from  <R,  4  R  )  data  near  r  =  r  (a) ,  and 

l  +  -  m 

0  ^ 

**2  from  (R+  -  R_)  data  near  r  =  rm(cr)  and  r  =  r  (a)  . 

/V  A 

f,  f  are  given  directly  by  (62) ^  and  then  F,  F  by  (62) 2,  and 

solution  over  a  maximum  (r,a)  domain  allows  complete  reproduction 

of  constant  stress  response  by  the  differential  law. 

In  the  reduced  model  (37),  (38),  the  argument  d2  of  F 

at  a  given  constant  a  is  duplicated  in  the  intervals 

r  4  r  4  r*  and  r*  <  r  <  r  where 
m  e 

r*  *  [?(rm  +  re>]'-  (64> 

Thus  the  relation  (62) 2  can  be  used  in  only  one  of  these  intervals. 

We  apply  (62)  on  the  intervals  r  £  r  «  r*,  r  £  r  *  r  ,  so 

m  e  o 

that  F  is  fully  determined,  and  (23) 2  on  r*  ^  r  $  r&  to 
complete  f.  Hence  primary  creep  is  matched  over  the  complete 
range  rm  <  r  <  rQ,  but  tertiary  creep  only  over  rm  r  $  r*, 
so  that  the  reduced  model  predicts  tertiary  creep  over  r*  $  r  <  r  , 
continuous  with  data  at  r*  and  rg.  The  important  primary-tertiary 
transition  is  matched.  This  correlation  has  determined 
f  ■  q (o)I{  d i)  and  F  =  q2(c)F(  d2)  on  0  <  <  *q(o)  “  r*(°> 

and  0  <  Cl2  <  ^Jr2(o)  “  rm^°^J  '  As  the  constant  o  is  varied 
there  will  be  a  common  range  of  and  of  d^  over  which  f 

and  F  are  already  determined,  and  hence  f  and  F  within  the 


2 

factors  q(o),  q  (o)  respectively.  But  independent 
response  R+(r,o),  R_(r,o)  at  a  different  a  will  require 
f(r ,o)  satisfying  (62)  with  r  dependence  in  general  not 
compatible  with  the  above  functions.  Thus  the  restricted  r 
dependence  of  (38)  allows  the  detailed  time  response  at  one 
constant  o  only  to  be  described,  but  does  incorporate  the 
properties  rm(o),  re(o),  and  rQ(a),  for  all  a  .  We 
suppose  rm(a),  r e ( a ) ,  rQ(a)  all  increase  with  o  ,  and  also 
that  rQ  -  rm  and  r*  -  rm  increase  so  that  the  ranges  of  both 
0. ^  and  CI2  increase  as  0  increases.  Hence  we  correlate 
the  time  response  for  the  maximum  a  ,  -  l  say.  Setting 
q(Z)  =  1,  ?(  G^)  and  P(  Q^)  are  now  fully  determined,  and 

the  factor  q(o)  allows  one  more  feature  of  the  constant 
o(<E)  curves  to  be  reproduced. 

For  illustrations  we  prescribe  in  turn  tm(cx),  related 
to  q(o)  by  (40),  and  em(o)  =  0*01  where 

e  Wo) 

cm  “  e  -  1,  cm(o)  =  r (t)dt .  (65) 

'  o 

By  (23),  (38)  r_  =  -  q(o)  { T (  ax )  +  p2(0tl)  +  which 

has  a  solution  for  each  constant  0  of  the  form 

r_(t)  *  e[q(o)t,o3  where  0(0, 0)  «  r0(o)  and  0(t*,o)  is 

determined  by  1  and  F,  independent  of  q(a).  By  (40), 

*  _  _ 
tm  (0)  =  q(o)tm(o)  is  determined  by  f,  F  ,  independent  of 

q(o),  so  changing  variable  in  (65)  to  t*  =  q(o)t  gives 

tMo) 

m 

e  “  qToT  I  eQt*,oldt*  .  (66) 

Thus  q(o)  is  determined  by  prescribing  e  (0),  hence  e. 
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since  the  integral  is  independent  of  q(o). 

For  each  constant  r  the  given  stress  response  determines 
o+(t)  =  S+(r,a)  >  0  and  a_  (t)  =  S_(r,  )  0,  and  by  (42) 

0  i  o  $  oM:  g  =  S+  +  S_,  (g2  +  4F) ^  =  S+  -  S_  ,  (67) 


with  F  e  0  and  S_  =  0  for  O  4  o  4  av  (r  r  )  .  Both 
relations  (67)  cannot  be  satisfied  since  F  has  been  determined 
by  constant  stress  response.  This  is  the  extent  to  which 
constant  o  and  constant  r  response  reflect  a  common 
material  property  under  this  model,  g  can  be  calculated  by  either 
(67 ) ^  or  (67)2,  when  the  other  fails  in  general,  or,  perhaps 
better,  the  complete  initial  response  OM  (Fig.  2)  can  be 
reproduced  by  matching  o ,  from  (42) : 

g  +  (g2  +  4F) ^  -  2S+  .  (68) 


The  time  variation  o(t)  during  relaxation  ME  (Fig.  2)  is 

then  predicted,  but  tM(r),  oE(r)  match  given  data. 

In  the  reduced  model  (54),  G  =  F  and  q  are  already 

known,  leaving  glC^)  and  s(r)  free.  Setting  s  =  1  at  the 

maximum  strain-rate  rQ(E) ,  g (  Q,^)  can  be  calculated  by  (66) 

to  reproduce  the  initial  response  OM  at  rQ(E).  Then  for 

r  (E)  ^  r  $  r  (E),  s(r)  can  be  determined  to  match  prescribed 
mo 

times  tM(r),  for  which  (48)  becomes 


Vr)  * 


1 


°Mlr) 


2do  1 


o  q(o  ){  [s2(r)g2(  aj)  +  4F  (  a  ^  %  +  s(r)g(fl,£)} 


(69) 
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where 

ai  -  r2  -  r*(o'),  &2  =  t2  -  rm(0')][re(0,>  -  *2]>  <70> 

or  to  match  the  strain  at  maximum  stress  eM(r)  =  e[oM(r)]. 

The  integral  is  evaluated  for  a  sequence  of  s (r) ,  each  r, 
until  t„(r)  is  attained,  in  our  illustrations,  since  we 

M 

cannot  construct  constant  strain-rate  responses  which  correspond 
to  the  idealised  constant  stress  responses,  we  set  g  =  f 
and  consider  both  t„(r)  =  t  [Pufrf]  and  e„(r)  =  0*01, 
that  is  eM  =  £n(l*01). 

Numerical  evaluation  of  the  complete  functions 
f(r,a),  g(r,o),  F(r,o),  by  (62),  (63)  and  (67)  or  (68)  in 

the  dependence  domain  of  the  (r,o)  plane,  shown  in  Fig.  3, 
is  straightforward,  but  the  large  matrices  of  tabulated  values 
are  worthwhile  only  when  a  full  set  of  good  constant  stress  and 
strain-rate  data  is  available. 

Idealised  response  and  model  correlation 

To  demonstrate  the  construction  of  the  simplified  reduced 
model  (37),  (38),  (54),  and  the  degree  to  which  it  describes 
the  complete  constant  stress  response,  we  have  treated  two 
different  families  of  idealised  curves  of  the  shape  shown  in 
Fig.  1.  Each  has  the  form 

r  =  r  (o)  -  Be“bt  +  Ce‘ct,  (71) 

where  the  stress  dependent  positive  parameters  B,  b,  C,  c 
are  chosen  to  match  given  r^(o),  r  (o),  x  =  t./t„  where  t. 

O  iu  X  In  1 


is  the  time  to  inflexion  (r  =  O)  on  the  tertiary  creep  curve 

me,  and  t  (o)  or  e  (o).  In  each  family  the  Colbeck  and 
mm 

Evans  (1973)  polynomial  flow  law  near  melting  is  adopted  for 
the  minimum  strain-rate: 

r  =  0*21o  +  0* 14o3  +  0*055o5,  (72) 

m 

5  -2 

where  o  is  measured  in  bars  (10  Mm  )  and  strain-rate  in 
years"1. 

In  family  1 : 

_  1 

re  *  rm(1  +  °h)t  ro  =  rm(1  +  2°H)t  fcm  =  tl°  T  "  1>8'  {73) 

where  t^  is  value  of  tm  at  o  =  1  bar.  To  be  specific  we 

adopt  t^  =  15  hours,  within  the  range  indicated  by  Colbeck 

and  Evans  (1973).  Choosing  the  maximum  stress  Z  ■  2  bars,  with 

q(E)  «■  1,  the  relations  (62)  determine  f(Cfc^),  F(  &2) ,  and 

q(o)  for  a  <  Z  is  determined  by  the  t  relation  (40).  Two 

m 

constant  strain-rate  responses  are  modelled  by  setting 
g(  dj^)  »  1(  (X^)  and  determining  s  (r)  =  s^fr)  and 
s (r)  »  s^fr)  in  the  g  decomposition  (54)  from  (69)  with 

tM(r)  *  tmC’(,M{r)l  "  b*  (2)'  and  tM(r)  *  Ft>M{r)]/r' 

respectively.  The  functions  7(  Q^),  F(&2),  X/q(a), 

l/s1(r),  l/s2(r),  are  shown  in  Fig.  4.  A  selection  of  constant 
stress  responses  r(t)  at  o  «  Oy  given  by  the  idealised 
functions  (72),  are  shown  by  the  solid  lines  in  Fig.  5,  and  the 
corresponding  responses  predicted  by  the  differential  operator 
model  (22)  are  shown  by  the  dashed  lines.  By  construction  the 
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o  »  E  =  2  bars  curves  match  identically  expect  on  the  final 

stage  r*  <  r  <  r  of  the  tertiary  creep,  and  r  ,  r  ,  r  ,  t 

e  oemm 

match  for  all  a  .  The  agreement  between  the  reduced  model 

and  the  parent  idealised  response  is  excellent.  Figure  6  shows 

the  corresponding  strain  responses  e(t)  calculated  from 

the  given  strain-rates  and  predicted  strain-rates.  Stress-strain 

curves  at  constant  strain- rates  r  =  r  (c.)  are  shown  in 

m  3 

Fig.  7  for  s  =  s  ( - )  and  s  =  s  2  ( - ).  Note  that 

a  differential  operator  law  which  predicts  primary  creep  at 
constant  stress  with  monotonic  r  still  allows  a  primary 
stress-strain  (time)  response  at  constant  r  with  an  inflexion 
point.  The  assumption  g  =  f  is  of  course  artificial,  and 
matching  with  given  constant  strain-rate  response  may  yield  a 
different  picture. 

In  family  2: 


re  ’  rmP-  +  0(o)l'  ro  *  rm&  +  ke<°0*  “  V"1'1 


9 


1*5, 

(74) 


and  0(a)  is  determined  by  setting  em  =  0*01  in  (65).  Again 
7  and  F  are  determined  by  matching  the  a  =  E  *  2  bars  response 
except  for  tertiary  creep  on  r*  <  r  <  rg.  First  we  determine 
q(o)  ■  q1(o)  by  matching  the  prescribed  (a  >  for  a  <  l 
through  the  relation  (40),  and  find  as  expected  that  the 
independent  property  em(°)  "  £<,(<*)»  calculated  from  the  predicted 
r(t),  is  not  0*01.  The  comparison  of  prescribed  ( — )  and 
predicted  ( - )  r(t)  curves  at  selected  stresses  a  »  a.  is 
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shown  in  Fig.  8,  and  the  comparison  of  strains  e (t)  determined 

and 

by  the  prescribed/_predicted  strain-rates  is  shown  in  Fig.  9, 
highlighting  the  differences  in  given  em  and  predicted 
em  *  ec  .  Next  we  determine  q  =  q2(o)  for  the  predicted 
strain  cm (a )  at  minimum  strain  rate,  given  by  (65),  (66), 
to  be  0*01,  and  show  the  corresponding  r(t)  and  c (t) 
curves  (....)  in  Figs  8  and  9.  Now  the  predicted  tm(o^) 
do  not  match  the  prescribed  t m (a  ^ ) .  The  primary  creep  prediction 
appears  better  for  the  choice  q  =  q2,  but  the  overall  match 
looks  better  with  q  =  q^.  Matching  of  the  reduced  model 
prediction  with  the  given  response  is  clearly  less  satisfactory 
for  family  2  than  family  1.  Similar  discrepancies  have  been 
obtained  with  t  =  1*6  and  k  =  5  in  (74).  A  distinct 
difference  between  the  two  families  is  that  6  ■  o*5  ■+■  0 
as  a  0  in  family  1,  while  6(o)  •*>  eQ  >  0  as  a  *  0  in 
family  2,  and  while  the  response  functions  7,  F,  s  are 
similar  for  the  two  families,  illustrated  in  Fig.  4,  1/q  ■+  0 

as  a  -*■  0  in  family  1,  but  not  in  family  2. 

Constant  strain-rate  responses  at  r  =  rm(Oj)  are 
constructed  by  determining  (i)  s  =  s^fr)  for 
tM(r)  -  t^r”1  (r)J,  (ii)  s  =  s2(r)  for  tM(r)  =  eCoM^ril/r, 
(iii)  s  *  Sj(r)  for  =  0*0l/r,  all  with  q  =  q^fo). 

These  are  illustrated  by  the  stress-strain  curves  shown 

( — — ) ,  (-—-),  ( . )  respectively  in  Fig.  10.  The  significant 

peaks  as  r  decreases,  that  is,  large  °M/0£*  arise  because 
r#/rw  •*  1  as  a  -*•  0  in  family  1,  but  approaches  11*45  in 
family  2. 
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Section  2.  Viscoelastic  solid  laws  of  differential  type 
Introduction 

The  viscoelastic  fluid  law  described  in  section  1  is 
necessarily  isotropic  in  all  configurations,  and  to  describe 
anisotropy,  whether  in  an  initial  reference  configuration  or 
in  a  subsequent  configuration  reached  by  loading,  a  solid 
law  is  required  with  explicit  dependence  on  a  reference 
configuration.  No  simple  non-linear  viscoelastic  solid 
laws  have  been  constructed  or  analysed  in  other  fields  of 
rheology,  so  we  adopt  a  heuristic  construction  by  analogy 
with  the  successful  fluid  law  of  differential  type. 

Eventually  we  hope  to  construct  an  integral  operator  law 
which  allows  more  general  dependence  on  deformation  history, 
and  which  should  be  more  satisfactory  for  numerical 
computation.  A  preliminary  examination  of  an  integral  law 
linear  in  finite  strain  shows  this  cannot  describe  important 
features  of  ice  response,  so  a  new  approach  must  be  investigated. 

Mellor  (1980)  notes  that  the  typical  constant  stress  and 
constant  strain-rate  responses  shown  in  Figs  1  and  2  are  often 
constant  load  and  constant  displacement- rate  responses, 
at  least  for  small  strains  achieved  in  laboratory  tests. 

This  is  an  immediate  reference  to  the  initial  configuration. 

For  uni-axial  compressive  stress  a  >  0  there  is  an  axial 
stretch  A  <  1  (contraction)  and  equal  lateral  stretches 
A-*  for  incompressible  material.  The  deformation  gradient 
and  associated  Cauchy-Green  strain  tensors  are 
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’  A 

0 

0 

A2  0  O' 

F  = 

0 

0 

»  B  =  C  = 

0  A_1  0 

0 

0 

X-11  ] 

0  0  A-1 

with  strain-invariants  det  B  5  1  and 
Kx  =  trB  -  X2  +  2A_1 ,  K2  =  i{(trB)2  -  trB2}  =  2A  +  A~2.  (2.2) 

The  axial  compressive  (engineering)  strain,  contraction  per 
unit  initial  length,  is 

e  =  1  -  A  ,  (2.3) 

and  the  compressive  strain-rate  relative  to  the  current 
configuration  is 

r  .  -A_1 A  =  (1  -  e)-1e  .  (2.4) 

Thus  constant  displacement-rate,  which  implies  constant  A  and 

e  ,  implies  constant  r  A  .  In  particular,  an  asymptotic 

strain-rate  re  >  0  implies  that  A  -*•  0  as  A-*-o,  e-*-l 

as  e  -*•  O.  With  small  instantaneous  elastic  strain  e  when 

e 

a  load  is  applied  and  maintained,  the  typical  monotonic  e  -  t 
curve  is  shown  in  Fig.  11,  indicating  a  1:1  relation  so  that  the 
engineering  strain-rate  e  may  be  expressed  as  a  function  of  e 
as  shown  in  Fig.  12.  The  dashed  sections  indicate  the  assumed 
finite  strain  response  consistent  with  the  above  asymptotic 
behaviour,  but  we  are  mainly  concerned  with  a  model  which 
describes  the  small  strain  response,  before  a  maximum  strain-rate 
and  subsequent  strain-rate  decrease  are  reached.  Mellor  suggests 


that  em#  corresponding  to  minimum  strain-rate 

em  =  (1  -  em)rm»  is  approximately  0*01  over  a  wide  range 

of  applied  stress. 

Force  per  unit  initial  (reference)  area  is  measured  by 
the  nominal  stress  o  =  F  *0  (det  F  =  1 ) ,  so  that 

constant  load  is  equivalent  to  constant  a  .  In  uni-axial 
compressive  stress,  with  =  -0, 

—  —  1  •  —  •  _ 

a  =  X  a  ,  a  =  Xo+  Aa  .  (2.5) 

At  constant  displacement  rate,  constant  e  ,  the  typical 
stress-strain  curve  is  shown  in  Fig.  13,  with  the  dashed 
section  indicating  response  beyond  small  strain 

observations  in  the  laboratory.  Stress  increases  to  a 
maximum  (?M  at  strain  eM  which  is  also  expected  to  be 
approximately  0*01  for  a  wide  range  of  e  ,  then  relaxes, 
but  as  compression  increases  there  must  be  a  minimum  load 
followed  by  unbounded  increase. 

The  most  simple  stress  dependence  on  a  reference 
configuration  is  through  a  frame  indifferent  tensor  function 
of  deformation  gradient  F  .  Then  in  uni-axial  stress,  o 
depends  explicitly  on  e  .  In  the  fluid  model  response  to 
constant  stress,  (18)  of  section  1,  it  was  necessary  to 
incorporate  a  strain-acceleration  term  r  to  exhibit  the 
time  variation  r(t)  shown  in  Fig.  1.  Here,  however, 
dependence  on  strain-rate  D  ,  and  hence  on  e  ,  is 
sufficient,  since  an  e,  e  relation  allows  the  strain-rate 
variation  shown  in  Fig.  12.  The  1:1  e  -  t  relation  shown 


33 


in  Fig.  11  replaces  the  time  variable  by  strain,  but  the 
differential  fluid  model  has  no  explicit  dependence  on  strain 
and  {18}  gives  constant  r  at  constant  o  if  the  r  term 
is  absent.  Again  we  require  a  stress-rate  term  to  exhibit 
a  time  variation  a{t),  and  hence  o(e)  variation 
shown  in  Fig.  13  at  constant  displacement  rate.  By  analogy 
with  (12)  in  section  (1),  an  appropriate  frame-indifferent 
differential  incompressible  solid  law  is 

S  +  *3fs(1)  -  |tr(S  D)l]  = 

[F  f  ( C )  FT  -  ^tr(C  f)l]  +  +  4>2[d2  “  fI2^  '  (2.6) 

where  f  is  an  arbitrary  symmetric  tensor  function  of 
symmetric  tensor  argument,  and  <J 4>2  are  functions 
of  any  invariants.  We  have  normalised  with  unit  stress 
deviator  coefficient,  and  assumed  linearity  in  the  stress 
tensor  rate 

If  the  solid  is  isotropic  in  the  reference  configuration 

then 

F  f  (C)  FT  -  ^tr(C  f)  =  [B  -  ^Kjl]  +  u2[B2  -  -j(K3  -  2K2)l],(2.7) 

where  u>2  are  functions  of  the  invariants  K^,  K2,  and 

possibly  of  J2»  J^,  and  the  rate  invariants.  Explicit 
dependence  on  e  in  uni-axial  response  can  be  retained  by 
including  dependence  of  <t>^,  <)>2  on  K]_  and  K2 
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*1  =  a,2  =  0  ;  that  is,  with  no  dependence  on  the  deformation 
orientation.  The  resulting  response  will  depend  on  the 
amount  of  deformation  from  a  reference  configuration,  but  can 
exhibit  no  anisotropy.  We  will  adopt  the  isotropic  form  (2.7) 
with  (2.6)  to  show  compatibility  with  uni-axial  response,  and 
we  plan  to  investigate  a  particular  form  to  demonstrate 
anisotropy  in  configurations  reached  by  simple  loading  paths. 
Initial  anisotropy  in  the  reference  configuration  arising 
from  the  ice  formation  will  require  a  different  shape  of  tensor 
function  f (C) . 


Uni-axial  response  and  small  strain  approximation 

In  uni-axial  stress  described  by  the  expressions 
(13)  -  (16)  of  section  1,  and  (2.1)  -  (2.5),  each  principal 
component  of  (2.6),  (2.7)  gives 


Xo  +  t^QXo  + 


2Xo] 


3, 

2*1 


(2.8) 


At  constant  load,  a  =  constant,  this  is  a  first  order 
differential  equation  for  A ( t )  subject  to  an  initial  condition 
X  (0)  **  1  -  eg(o).  From  Fig.  11  we  see  that  X  (t)  is  monotonic 
with  X  a  single-valued  function  of  X  .  This  suggests  we 
consider  a  simplified  form  of  (2.8)  with  a  absent,  by 
setting 


2 


0 


9 


(2.9) 


though  (2.9)  is  not  a  necessary  condition  for  one  real  root 
X  on  0  <  X  <  1.  At  constant  displacement  rate,  X  =  constant, 
we  have  a  first  order  differential  equation  for  o(t)  with 
X  =  At  and  a (0)  =  0.  The  necessary  a  term  to  exhibit 
the  time  variation  of  Fig.  13  implies  fj  ?  0,  though 
stress-rate  dependence  could  instead  be  incorporated  by 

9  9 

dependence  of  other  coefficients  on  J2»  J^,  with  i|<3  =  0. 

Figure  14  shows  the  variation  of  o  with  e  corresponding 
to  the  response  in  Fig.  13. 

A  small  strain  approximation  neglecting  |e|  compared 
to  unity  gives  lead  order  expressions  for  uni-axial  response 

1  -  X  =  e  ,  -X  =  e~r,  K.  -  3  =  3e2,  K2  -  3  -  3e2,  (2.10) 

and  ,w2  can  be  regarded  as  functions  of  o,  e  if 

rate-dependence  is  restricted  to  the  terms  and  D. 

The  observed  non-linear  response  at  constant  a  implies 
that  the  coefficients  vary  significantly  with 

small  changes  of  e,  so  cannot  be  approximated  by  their 
(constant)  values  at  =  K2  =  3.  To  compare  the  magnitudes 
of  various  terms  in  (2.8)  it  is  helpful  to  introduce  a 
normalised  strain  and  normalised  time  by 

e  -  Ce,  t  =  t/tQ  (2.11) 

where  C_1  is  a  maximum  strain  (*05  say)  and 

|o|,  —  and  —I  are  order  unity  in  magnitude  at  their 

•dt1  'dt1 

maxima.  We  suppose  the  stress  unit  is  one  bar  and  consider 
a  range  |o)  <  10  bars.  Now  (2.8),  (2.9)  reduce  to 
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~  ~  -  C_1fe(3u)1  +  6u>2)  +  ^(2 t^o  +  l^f]  =  0, 

o  dt  o  dt  -1 


o  + 


(2.12) 


and  immediately  |C  ^  a  de/dtj  <<  da/dt  ,  but  to  retain 

the  required  dependence  on  S^,  D,  and  B  ,  it  is  necessary 

to  define  order  unity  coefficients 


__  4-  _  3$.  _  3o).  +  6a)- 

=  ?1=2TT  '  =  =  — c— 

o 


(2.13) 


so  that  (2.12)  becomes 


—  t  do - 

o  +  i|i  —  -  id  e 

J  dt 


-4,^  =  0 


1  dt 


(2.14) 


Uni-axial  response  cannot  distinguish  the  w.  ,(d2  terms,  so 
we  have  three  independent  coefficients  w(o,e),  ^(o^e), 

43 (o,e) . 

The  differential  relation  (2.14)  does  not  apply  directly 
to  a  stress  jump  Aq  ,  but  we  can  deduce  the  resulting  strain 
jump  Ae  if  we  assume  that  (2.14)  applies  in  the  limit  of 
a  smooth  rapid  stress  change  from  oQ  to  oQ  +  Aq  and  also 


that 


e  =  ec  +  h (o)  ,  h(oQ)  =  0  , 


(2.15) 


for  any  positive  AQ  ,  accompanied  by  a  strain  change  eQ 

to  i?  +  A  where  A^  =  h  (o  +  A  ) .  Let  the  change  take  place 
o  e  e  o  o 

in  the  time  interval  (tQ,  t  +  6 1 ) ,  integrate  (2.14)  over  the 

time  interval,  and  let  St  -*■  0.  Since  uJ  is  bounded  we  obtain 

o+A  c  +A_ 

o  o  o  a 

j_  43  |jj,eo  *  h(a)Jdo  »  j_  h' (a)^^,^  +  h(a)jda,  (2. 


16) 
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for  all  Aa  ,  and  hence 

[a,eQ  +  h(o)]h'(o)  -  ?3£o,eo  +  h(a)J  =  0,  h laQ)  -  0.  (2.17) 


Thus,  given  4> , ,  ip  ,  (2.17)  is  a  first  order  differential 
i  3  for 

equation  with  initial  condition/^ h (a) ,  and  the  solution  determines 

Ae  =  h  (cjo  +  A^)  for  all  Aq  .  The  strain  jump  depends  on  the 

stress  jump  and  the  starting  stress.  In  particular,  if  the 

jump  A  is  applied  at  time  t  =  0  from  zero  stress  and  strain, 

Ae  is  given  by  setting  eQ  =  oQ  *  0  in  (2.17),  and  determines 

the  elastic  strain  e  (A)  =  C  ^A  (A).  However,  the  uni-axial 

e  e 

response  for  all  A  determines  A0(A)  =  h(A),  and  so  in 
principle  determines  a  relation  between  and  ij>3  in  the 

strain  range  e  =*  O  -*•  C  e  not  covered  by  the  smooth  response. 


Constant  load  and  constant  displacement  rate. 


For  constant 


in  t  >  0,  (2.14)  becomes 


i  -  +  we  -  a  =  0  .  (2.18) 

1  dt 

Given  a  family  of  normalised  curves  corresponding  to  Fig.  12, 


—  «  F(o,e)  ,  (2.19) 

dt 

then 

o  -  we  =  $^F  ,  e  i  eg(o) 


(2.20) 
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determines  one  relation  between  u  and  4>^.  For  constant 
displacement  rate  defined  by  de/dt  =  w,  so  e  =  wt  ,  (2.14) 
becomes 

rJifT 

—  +  a  -  a)  e  -  $.w  =  0  .  (2.21) 

*  dt  L 

Given  families  of  normalised  curves  corresponding  to  Figs  13  and 
14/ 

a  =  g(w,e),  —  =  w  =  G(w,e),  e  0  ,  (2.22) 

dt  de 

then 

^3 (w#e)  ,ej  G  (w,e)  =  [g  (w,e)  ,e]{-F  (g(w,e)  ,ej  +  w}  (2.23) 

determines  in  terms  of  .  Since  F  >  O,  w  >  O,  and 

G  =  O  at  e  -  eM(w),  and  assuming  ^  (o  ,e)  has  no  zero,  then 

F[g{w,eM(w)  },  eM(w)]  =  w  .  (2.24) 

Now  (2.20),  (2.23)  provide  two  relations  for  three 
functions  ijL,  w,  all  of  which  are  required  to  describe  the 

responses  in  Figs  11  and  13  and  incorporate  dependence  on  the 
reference  configuration  through  B  .  Clearly  w  =  0,  o  =  ^F  , 
solve  (2.20),  but  induced  anisotropy  is  no  longer  possible. 

In  this  case,  however,  ^  and  lj>3  are  determined  by  F  and 
G,  and,  alternatively,  complete  responses  F  and  G  can  in 
principle  be  matched  by  suitable  choice  of  and  When 

u  /  0  there  is  additional  flexibility. 


While  some  values  of  g  =  o  are  repeated  on  the  two 
sides  of  e  =  eM(w),  illustrated  in  Fig.  13,  the  arguments 
(g,e)  corresponding  to  the  same  g  on  a  constant  w  curve 
are  distinct  points  of  the  (o,e)  domain.  Further,  assuming 
that  3g/3w|—  >  0,  the  family  of  curves  are  non-intersecting, 
so  each  constant  w  relation  (2.23)  between  ^  and  4^ 
refers  to  distinct  sets  of  points  (o,e) .  This  situation 
contrasts  strongly  with  the  minimal  fluid  model  in  section  1 
which  was  not  sufficiently  general  to  match  independent 
constant  stress  and  constant  strain- rate  data. 

A  restricted  model  =  ^(e),  u  =  io(e)  satisfies 

(2.20)  provided  that 

F  =*  ofx(i)  +  fQ (e) ,  *x  =  l/fr  w  =  -^1fQ/e  ,  (2.25) 

which  requires  the  strain-rate  at  fixed  strain  to  change 
linearly  with  the  constant  stress  o  .  A  positive  f1(e) 
would  imply  a  linearly  increasing  strain-rate  which  we  expect 
to  be  qualitatively  correct.  With  the  restrictions  (2.25), 

(2.23)  becomes 

?3j9<w»e)  ,<fjG(w,e)f (e)  =w  -  fQ  (e)  -  f  x  (e)g  (w,e) ,  (2.26) 

which  requires  that  ipj  depends  on  both  o  and  e  in  general. 
Writing  ,  (2.23)  becomes 

-  e(s(w,5),5]u  -  sksnuiLl) )  (2.27) 

3e 

which  interprets  E(o,e)  as  a  Young's  modulus  in  the  limit 
of  inifinite  strain-rate.  A  restriction  E  *  constant  Implies 
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a  relation  between  the  supposed  independent  functions 
F (o,e),  g (w,e) >  that  is  from  (2.22),  (2.27), 


3G|  _  _  E  f  G  3F  +3F.  ) 

3e  w  l w  3o  e  3e  o  ' 


(2.26) 


If  e  *  e^(a)  corresponds  to  the  second  inflexion  point  in 
Fig.  11,  or  local  maximum  e  in  Fig.  12,  and  e  =  e2  (w) 
corresponds  to  the  inflexion  point  in  Fig.  13,  or  local  minimum 

e 

o  in  Fig.  14,  then 


3F 

3e 


<  0 
>  0 
<  o 


0  <  e  <  e_ 


m 


e„<  e  <  e. 
m  1 


eL  <  e  , 


3G 

3e 


w 


<  0 
>  0 


0  <  e  <  e2 


e2 


<  e 


(2.29) 


Adopting  the  anticipated  inequalities  3F/3oi_  >  0  and 

I e 

«2  >  em  '  then  the  inequalities  (2.29)  are  certainly  consistent 
with  (2.28)  for  em  ^  e  ^  m*n  (®1'®2)'  but  no  ^irm  conclusion 
can  be  drawn  in  the  other  ranges  of  e  .  Thus  E  =  constant 
may,  or  may  not,  be  an  acceptable  model.  Given  independent 
F,g,  (2.27)  determines  E(a,e)  ,  and  (2.20)  relates  w  to 

\  or  *3  • 

The  next  stage  of  our  investigation  must  explore  various 
simplications  to  find  u,  satisfying  (2.20),  (2.23), 

for  F,g,G  compatible  with  expected  qualitative  response. 

Once  such  a  model  is  determined  we  wish  to  demonstrate  the 
anisotropy  induced  by  simple  loading  histories. 
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Figure  Captions 

1.  Strain-rate  r(t)  at  constant  stress  o  . 

2.  Stress  o(t)  at  constant  strain-rate  r. 

3.  Dependence  domains  of  response  functions. 

4.  Functions  f(CL^),  F (ft2),  l/q(CT),  1/s^r),  l/s2(r)  for 
family  1. 

5.  Strain-rates  r(t)  at  constant  o  =  ck  =  2 (-0* 25)0* 5  bars, 

for  family  1  ( - ) ,  predicted  by  reduced  model  ( - ) . 

r  decreases  at  fixed  t  as  decreases. 

6.  Strains  e(t)  corresponding  to  strain-rates  shown  in 
Fig.  5. 

7.  Stress-strain  curves  at  constant  r  =  rm(o..)  predicted 

by  reduced  model  for  family  1  with  s  =  s^  ( - )  and 

s  »  s2  ( - ) . 

8.  Strain-rates  r(t)  at  constant  a  =  =  2 (-0*25)0*75  bars, 

for  family  2  ( - ) ,  predicted  by  reduced  model  with 

q  ■  ( - )  and  with  q  =  q2  ( . )• 

9.  Strains  e(t)  corresponding  to  strain-rates  shown  in  Fig.  8. 

10.  Stress-strain  curves  at  constant  r  =  rm(o..)  predicted  by 

reduced  model  for  family  2  with  q  =  q^  and  s  *  s^( - ), 

s  ■  S2  ( - )  ,  s  =  S3  ( . )  . 

11.  Compressive  strain  response  at  constant  load. 

12.  Strain-rate  v.  strain  at  constant  load. 

13.  Stress-strain  curve  at  constant  displacement-rate. 

14.  Stress-rate  v.  strain  at  constant  displacement-rate. 


rate  r(t)  at  constant  stress 


Functions 


^4l>»  F(CL2),  1/q (o)  ,  1/s^r),  l/i 


(M)  09 


Strains  e(t)  corresponding  to  strain-rates  shown  in  Fig. 


Fig.  7.  Stress-strain  curves  at  constant  r  =  r  (o.)  predicted  by  reduced  model  for 


Compressive  strain  response  at  constant  load 


Pig.  12.  Strain-rate  v.  strain 


Fig.  14.  Stress-rate  v.  strain  at  constant  displac* 


